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We study the nucleon and diquarks in a non-local Nambu-Jona-Lasinio model. For certain pa- 
rameters the model exhibits quark confinement, in the form of a propagator without real poles. 
After truncation of the two-body channels to the scalar and axial-vector diquarks, a relativistic 
Faddeev equation for nucleon bound states is solved in the covariant diquark-quark picture. The 
dependence of the nucleon mass on diquark masses is studied in detail. We find parameters that 
. - . lead to a simultaneous reasonable description of pions and nucleons. Both the diquarks contribute 

\l ' attractively to the nucleon mass. Axial- vector diquark correlations are seen to be important, espe- 

I cially in the confining phase of the model. We study the possible implications of quark confinement 

^ . for the description of the diquarks and the nucleon. In particular, we find that it leads to a more 

' compact nucleon. 
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,— I ■ I. INTRODUCTION 

> : 

' The NJL model is a successful phenomenological field theory inspired by QCD T|. The model is constructed to obey 
, the basic symmetries of QCD in the quark sector, but unlike the case of low-energy QCD, quarks are not confined. 
QQ ' The basic ingredient of the model is a zero-range interaction containing four fermion fields. This means that the 
, model is not renormalizable. Therefore at one-loop level an ultraviolet cut-off supplemented with a regularization 
■rj" ■ method is required from the outsets. The value of the cut-off can be related to the scale of physical processes not 
\ included in the model, and thus determines its range of validity. Consequently, processes involving a large momentum 
■ transfer can not be described by the model. At higher orders in the loop expansion, which are necessary for calculating 
rS i' mesonic (baryonic) fluctuations 2, 3], one needs extra cut-off parameters. It is hard to determine these parameters 
JL from independent physics, and thus to build a viable phenomenology. A similar problem appears in the diquark-quark 

picture of baryons where an additional cut-off parameter is required to regularise the diquark-quark loops |[3| . 
(~| Another drawback of the model is the absence of confinement, which makes it questionable for the description of 

I!" few-quark states and for quark matter. If energetically allowed, the mesons of the model can decay into free quark- 
. !^ antiquark pairs, and the presence of unphysical channels is another limitation on the applicability of NJL model. At 
the same time, it is also known that the NJL model exhibits a zero-temperature phase transition at unrealistically 
5-H ' low baryon density This problem is caused by the formation of unphysical coloured diquark states. These may be 
explicitly excluded at zero density by a projection onto the physical channels, but dominate the behaviour at finite 
density. The model is not able to describe nuclear matter, even in the low-density regime . 

We do not know how to implement colour confinement in the model and, anyway, the exact confining mechanism 
of QCD is still unknown. In the context of an effective quark theory, a slightly different mechanism of "quark 
confinement" can be described by a quark propagator which vanishes due to infra-red singularities or which does 
not produce any poles corresponding to asymptotic quark states 0, 0- Another realisation of quark confinement 
can be found in Ref. B. I t has been shown that a non-local covariant extension of the NJL model inspired by the 
instanton liquid model 'lO'| can lead to quark confinement for acceptable values of the parameters This model 
has previously been applied to mesons [llL \i% .13.] and in this paper it is applied to baryons based on the relativistic 
Faddeev approach. 

The quark propagator in the model has no real pole and consequently quarks do not appear as asymptotic states. 
Instead the quark propagator has pairs of complex poles. This phenomenon was also noticed in Schwinger-Dyson 
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equation studies in QED and QCD IMlialig- 

One can simply accept the appearance of these poles as an artifact 
of the naive truncation scheme involved. However, it has been recently suggested that it might be a genuine feature 
of the full theory, and be connected with the underlying confinement mechanism . For example, it has been 

shown by Maris that if one removes the confining potential in QED in 2+lD the mass singularities are located almost 
on the time axis, and if there is a confining potential, the mass-like singularities move from the time axis to complex 
momenta 16]. In this paper, we study this kind of confinement from another viewpoint. We show that when we have 
quark confinement in the non-local NJL model, the baryons become more compact, compared to a situation where 
we have only real poles for quark propagator. 

There are several other advantages of the non-local version of the model over the local NJL model: the dynamical 
quark mass is momentum-dependent, as also found in lattice simulations of QCD pJl]. There various methods are 
available for construction of a conserved current in the presence of non-local interactions [T^ . In gen eral, one can 
preserve the gauge invariance and anomalies by introducing additional non-local terms in the currents [iSj . A Noether- 
like method of construction for these non-local pieces for the non-local NJL model was developed in Ref. The 
regulator makes the theory finite to all orders in the loop expansion and leads to small next-to-leading order corrections 
As a result, the non-local version of the NJL model should have more predictive power. 

We use a separable non-local interactions, similar to that of the instanton-liquid model |10l.|l9|. This considerably 
simplifies the calculation. Other approaches also give non-locality but in different forms 0, . Non- locality also 
emerges naturally in the Schwinger-Dyson resummation and in various types of gluonic field configuration within 
the QCD vacuum, see for an example Ref. [2H. 

Considerable work has been done on these nonlocal NJL models including applications to the mesonic sector 
[iTl I13I I , phase transitions at finite temperature and densities , and the study of chiral solitons j23| . 

In this paper we present our first results from a calculation of the relativistic Faddeev e quat ion for a non-local NJL 
model, based on the covariant diquark-quark picture of baryons [23. 125I |2^ WA . [281 12^ Isol |3l|j3|LSuch an approach 
has been extensively employed to study baryons in the local NJL model, see, e.g., Refs. [24 l25l l26l I271 128^ . We 
include both scalar and the axial-vector diquark correlations. We do not assume a special form for the interaction 
Lagrangian, but we rather treat the coupling in the diquark channels as free parameters and consider the range of 
coupling strengths which lead to a reasonable description of the nucleon. We construct diquark and nucleon solutions 
and study the possible implications of the quark confinement for the solutions. The dependence of the baryon masses 
and waves on the diquarks parameters is investigated and the role of diquarks in the nucleon solutions, for both the 
confining and the non-confining phase of the model is considered separately. The nucleon wave function is studied in 
details. Due to the separability of the non-local interaction, the Faddeev equations can be reduced to a set of effective 
Bethe-Salpeter equations. This makes it possible to adopt the numerical method developed for such problems in 
Refs. UllSllllll- 

This paper is organised as follows: In Sec. II the model is introduced. We also discuss the pionic sector of the 
model and fix the parameters. In Sec. Ill the diquark problem is solved and discussed. In Sec. IV the three-body 
problem based on diquark-quark picture is investigated. The numerical technique involved in solving the effective 
Bethe-Salpeter equation is given and the results for three-body sector are presented. Finally, a summary and outlook 
is given in Sec. V. 



II. A NON-LOCAL NJL MODEL 



We consider a non-local NJL model Lagrangian with SU{2)f x SU{3)c symmetry. 

C = ^ii^-m,)ij + Ci, (1) 

where rric is the current quark mass of the u and d quarks and Cj is a chirally invariant non-local interaction 
Lagrangian. Here we restrict the interaction terms to four-quark interaction vertices. 

There exist several versions of such non-local NJL models. Regardless of what version is chosen, by a Fierz 
transformation one can rewrite the interaction in either the quark-antiquark or quark-quark channels. We therefore 
use the interaction strengths in those channels as independent parameters. For simplicity we truncate the mesonic 
channels to the scalar (0+, T = 0) and pseudoscalar (0~, T = 1) ones. The quark-quark interaction is truncated to the 
scalar (0+, T = 0) and axial vector T = 1) colour 3 quark-quark channels (the colour 6 channels do not contribute 
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to the colourless three-quark state considered here) . We parametrise the relevant part of interaction Lagrangian as 

Cl = \g-^ia{x)ia{x) + gsJs{x)Js{x) + gaJa{x)Ja{x), 

x)^{xi)[ToCT2l3^]'f [x^), 

x)4,'^{x2)[C'^l^r2l3^]i^{xi), 

x)llj{xi)[-i^CTiT2l3^]i^^ {X^), 

x)^'^{x2)[C'^rr2nl3^]i^{x^), (2) 

where Ta = (1, *75t). The matrices (3^ = ■y/3/2A'^(A = 2, 5, 7) project onto the colour 3 channel with normalisation 
ir{j3^l3^ ) = iS^^ and the r^'s are flavour SU{2) matrices with tr(rirj) = 25.^. The object C = 17275 is the charge 
conjugation matrix. It is exactly this four-way separability of the non-local interaction that is also present in the 
instanton Hquid model [T^ . 

Since we do not restrict ourselves to a specific choice of underlying interaction, we shall treat the couplings gs, 
ga and iJtt as independent parameters. We assume g,r,s,a > 0, which leads to attraction in the given channels (and 
repulsion in the quark-antiquark colour octet and quark-quark colour antisextet channels). The coupling parameter 
is responsible for the pions and their isoscalar partner a. The coupling strengths gs and ga specify the behaviour 
in the scalar and axial-vector diquark channel, respectively. 

We define the Fourier transform of the form factor by 

^("^""^^^^ / (3) 

The dressed quark propagator S{k) is now constructed by means of a Schwinger-Dyson equation (SDE) in the rainbow- 
ladder approximation. Thus the dynamical constituent quark mass, arising from spontaneously broken chiral symme- 
try, is obtained in Hartree approximation as^ 

m{p) = + ig^f{p) j TT[S{k)]f{k), (4) 

where 



ja{x) = 


1 d'^xid'^X3F{x 


- X3)Fixi 


J4x) = 


j d^xid'^X3F{x 


- X3)F{xi 


Js{x) = 


j d'^X2d'^XiF{x 


- X4)F{X2 


Jaix) = 


j d'^xid'^x^Fix 


- X3)F{xi 


Ja{x) = 


1 d'^X2d'^X4F{x 


- X4:)F{X2 



S-\k)^l^~m{k). (5) 

One can simplify this equation by writing to(p) in the form 

to(p) = TOc + (r7i(0) - mc)/^(p). (6) 

The non- linear equation can then be solved iteratively for m(0). 

Following Ref. 'l^, we choose the form factor to be Gaussian in Euclidean space, /{pe) — exp(— p|./A^), where pe 
denotes the Euclidean four-momentum and A is a cutoff of the theory. This choice respects Poincare invariance and 
for certain values of the parameters it leads to quark, but not colour, confinement. For values of to(0) satisfying 

m{0)-mc 1 / (ymfTA2 + mc)2 

^y^+W^^c ^ 2 """"^ y 2A2 

the dressed quark propagator has no poles at real p^ in Minkowski space (p^ + m?{p'^) ^ 0). The propagator has 
many pairs of complex poles, both for confining and non-confining parameter sets. This is a feature of these models 




The symbol Tr denotes a trace over flavour, colour and Dirac indices and tr/3 denotes a trace over Dirac indices only. 





FIG. 1: A graphical representation of the Bethe-Salpeter equation for the qq T-matrix in RPA approximation. The sohd hnes 
denote the dressed quark propagators Eq. and shaded boxes denote meson propagators. 
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FIG. 2: One-pion-to-vacuum matrix element in RPA, contributing to the weak pion decay. The lines are as defined in Fig. 
The wavy line denotes a weak decay. 



and due care should be taken in handling such poles, which can not be associated with asymptotic states if the theory 
is to satisfy unitarity. One should note that the positions of these poles depend on the details of the chosen form 
factor and the cut-off, hence one may regard them as a pathology of the regularisation scheme. Since the choice of the 
cut-off is closely related to the truncation of the mesonic channels, (for example, if one allows mixing of channels, the 
cut-off and the positions of poles will change). Even though the confinement in this model has no direct connection to 
the special properties of the pion, there is an indirect connection through the determination of the parameters from 
the pionic properties. 

The quark-antiquark T-matrix in the pseudoscalar channel can be solved by using the Bethe-Salpeter equation in 
the random phase approximation (RPA), as shown in Fig.Q] see Ref. pj|. 



T{pi,P2,P3,P4) = /(Pl)/(P2)[«75 



1 + gixJixiq^) 
^5{pi +P2 -P3-P4), 



[«75Ti]/(p3)/(P4) 



(8) 



where 



J^(g ) = iTi 
= 6i 



(2^ 



.f\k)-imS{k)ronS{q + k)f{q + k), 



trz,[755(fc)755(fc + qW{k)f\q + fc), 



(9) 



where q denotes the total momentum of the quark-antiquark pair. The pion mass corresponds to the pole of 
T-matrix. One immediately finds that m.^- = if the current quark mass is zero, in accordance with Goldstone's 
theorem. The residue of the T-matrix at this pole has the form 



V''{jPl,P2) = i9^qq[^c ®t''® 75]/(Pl)/(P2), 



(10) 



where gT^qq is the pion-quark-antiquark coupling constant and is related to the corresponding loop integral by 



-2 



(11) 



The pion decay constant is obtained from the c oupling of the pion to the axial- vector current. Notice that due to 
the non-locality the axial- vector current is modified [ij, [13 and consequently the one-pion-to-vacuum matrix element 
gets the additional contribution shown in Fig. |5] This extra term is essential in order to maintain Gell-Mann-Oakes- 
Renner relation |llj and makes a significant contribution. The pion decay constant is given by 



ml 



^Tr[^75y(5(p-))75ra(5(p+))]/(p_)/(p, 



jTT[S{k)] 



d-^k 



Tr[y-(p_,p+)5(p_)75ra^(p+)] 



2ml J (27r)4 ' ^ " J {2itY 
X [f(k) {f{p+) + f{p^)) f{p+)f{p-)f{k) {f{k + q) + f{k q))], 



(12) 
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TABLE I: The parameters for the sets A and B, fitted to = 92.4 MeV and m,^ = 139.6 MeV. The resuhing values of the 
dynamical quark mass m(0) are also shown. 

Parameter set A set B 

m(0) (MeV) 297.9 351.6 

mo(0) (MeV) 250 300 

mc (MeV) 7.9 11.13 

A (MeV) 1046.8 847.8 

g^(GeV"^) 31.6 55.80 



TABLE IL The first two sets of poles of the quark propagator (in magnitude) in the Minkowski frame. 



set A 


set B 


±391 MeV 


±408 ± 238i MeV 


±675 MeV 


±1575 ± 307i MeV 



where VTr{p-,P+) is defined in Eq. 110(1 . with the notation p± — pi: ^q. 

The loop integrations in Eqs. (|9I12|) are evaluated in Euclidean space^. For the current model, the usual analytic 
continuation of amplitudes from Euclidean to Minkowski space can not be used. This is due to the fact that quark 
propagators of the model contain many poles at complex energies leading to opening of a threshold for decay of a 
meson into other unphysical states. Any theory of this type needs an alternative continuation prescription consistent 
with unitarity and macrocausality. Let us define a fictitious two-body threshold as twice the real part of the first pole 
of the dressed quark propagator m^. For a confining parameter set, each quark propagator has a pair of complex- 
conjugate poles. Above the two-body pseudo-threshold < — 4(m^)^, where q is the meson momentum, the first pair 
of complex poles of the quark propagator has a chance to cross the real axis. According to the Cutkosky prescription 
|33|. if one is to preserve the unitarity and the microcausality, the integration contour should be pinched at that 
point. In this way, one can ensure that there is no spurious quark-antiquark production threshold, for energies below 
the next pseudo-threshold, i.e. twice the real part of the second pole of the quark propagator. Note that it has been 
shown 34] that the removal of the quark-antiquark pseudo-threshold is closely related to the existence of complex 
poles in the form of complex-conjugate pairs. Since there is no unique analytical continuation method available for 
such problems, any method must be regarded as a part of the model assumptions |Tll ITsl Is^. Here, we follow the 
method used in Ref. |TH . 

Our model contains five parameters: the current quark mass rric, the cutoff (A), the coupling constants , and 
ga- We fix the first three to give a pion mass of TOtt = 139.6 MeV with decay constant — 92.4 MeV, while we take 
the value of the zero-momentum quark mass in the chiral limit mo(0) as an input. We analyse two sets of parameters, 
as indicated in Table ^ where set A is a non-confining parameter set, while set B leads to quark confinement (i.e., it 
satisfies the condition Eq. O). The position of the quark poles are given in Table ITTI The real part of the first pole of 
the dressed quark propagator can be considered in much the same as the quark mass in the ordinary NJL model. 
Since we do not believe in on-shell quarks or quark resonances, this is also a measure for a limit on the validity of the 
theory. The mass is larger than the constituent quark mass at zero momentum m(0), as can be seen in TableQ] As 
we will see appears as an important parameter in the diquark and nucleon solution, rather than the constituent 
quark mass. The same feature has been seen in the studies of the soliton in this model, where determines the 
stability of the soliton *23'|. The parameters gs and ga will be treated here as free parameters, which allows us to 
analyse baryon solutions in terms of a complete set of couplings. This is permissible as long as the interactions in 
Lagrangian are not fixed by some underlying theory via a Fierz transformation^. The coupling-constant dependence 
is expressed through the ratios = gsl g-K and Ta = ga/g-n- 

The quark condensate {iptp) = i Tr S'(O) is closely related to the gap equation, Eq. |@J. In the latter there appears an 



^ We work in Euclidean space with metric g'^" = 5'"' and a hermitian basis of Dirac matrices {7^ ,71/} = 25^,^ , with a standard transcription 

rules from Minkowski to Euclidean momentum space: fc" — > ik4, k^'^ — > —k^ 
^ Notice as well that the Hartree-Fock approximation is equivalent to the Hartree approximation with properly redefining coupling 

constants. Therefore, the Hartree approximation here is as good as the Hartree-Fock one, since the interaction terms are not fixed by a 

Fierz transformation. 
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extra form factor inside the loop integral. The quark condensate in the chiral limit is —(207 MeV)'^ and —(186 MeV)^ 
for sets A and B, respectively. These values fall within the limits extracted from QCD sum rules [s^ and lattice 
calculations jH^j having in mind that QCD condensate is a renormalised and scale-dependent quantity. In contrast 
to the local NJL model, here the dynamical quark mass Eq. ^ is momentum dependent and follows a trend similar 
to that estimated from lattice simulations [l3 • 



III. DIQUARK CHANNELS 



In the rainbow-ladder approximation the scalar quark-quark T-matrix can be calculated from a very similar diagram 
to that shown in Fig. ^ (the only change is that the anti-quark must be replaced by a quark with opposite momentum). 
It can be written as 

T{pi,P2,P3,Pi) = /(Pi)/(P2) [75^X2/3^] r(q)[C-i75r2/3^]/(P3)/(P4) 

X(5(pi +P2 -P3 -P4), (13) 

with 

r{ci) = TtHt^V (14) 
where q — pi + P2 ^ P3 + Pi is the total momentum of the quark-quark pair, and 

Js{q^) = ^Trl ^j\-k)[^,Cr2P^]S{-kf[C-'-t,T2f3^]S{q + k)f{q + k), 

= 6z y ^irD[-l,S{k)^,S{k + q)]f{k)f{q + k). (15) 

In the above equation the quark propagator S{k) is the solution of the rainbow SDE Eq. |(SJ|. The denominator of 
Eq. lO is the same as in the expression for the pion channel, Eq. JSJ, if gs = Qtt- One may thus conclude that at 
rs = 1 the diquark and pion are degenerate. This puts an upper limit to the choice of r^, since diquarks should not 
condense in vacuum. One can approximate T(g) by an effective diquark "exchange" between the external quarks, and 
parametrise T{q) near the pole as 

T{q) = 2igl,^^V'{q)D{q), D'^q) ^ - All, (16) 

where Mds is the scalar diquark mass, defined as the position of the pole of T{q). The strength of the on-shell coupling 
of scalar diquark to quarks, gdsqq is related to the polarisation Js by 

and V^{q) is the ratio between the exact T-matrix and on-shell (one-pole) approximation and describes the off-shell 
correction of the T-matrix around the diquark solutions [s^ • For the "on-shell approximation" we have (g) = 1 
and by definition on the mass shell q^ = one has V^{q)\g2^M2 ~ ^' 
There is no mixing between the axial-vector diquark and other channels, and so one can write the axial-vector 
diquark T-matrix in a similar form 

T(pi,p2,P3,P4) = /(pi)/(p2)[7MC^T,T2/3^]T^''(g)[C-l7.T2T,/3^]/(p3)/(p4) 

X(5(pi +P2 -P3 -P4), (18) 

with 

^'l)-'3'^'[l + g,jT(^q2)+l+g^jL^q2^\' (1^) 

where we decompose the axial polarisation tensor into longitudinal and transverse components: 

Jr(g') = ^Trl ^f{-k)[j^CnT2l3^]S{-kf[C-'rT2np^]S{q + k)f{q + k), 

= 6*1 -0^trn[rSik)rSik + q)]f{k)f{q + k) 

= Jl{q')i9'"' ' q'q" /q') + J^{q^)q''q'' /q^- (20) 
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FIG. 3: The denominator of the diquark T matrix for the longitudinal and transverse axial vector channel, for parameter set 
A at = 0.44. Note that there is no longitudinal pole. 
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FIG. 4: The scalar and axial- vector diquark mass as a function of Vs and r^, respectively, for both parameter sets. The dotted 
and the dash-dotted lines denote the quark-quark pseudo-threshold for set A and B, respectively. 



Wc find that the longitudinal channel does not produce a pole (see Fig.O, and thus the bound axial- vector diquark 
solution corresponds to a pole of the transverse T-matrix. The transverse component of T^'^^q) matrix can be 
parametrized as, 



T^%q) ^ 2^gi,^^V'^{q)D^%q), D^%q) = ^"^f/^' , (21) 



where V^lq) includes the off-shell contribution to the axial- vector T-matrix. The coupling constant gdaqq is related 
to the residue at the pole of the T-matrix, 

9daqq= ^\q-=Ml^- (22) 



A. Diquark Solution 



The loop integrations in Eqs. H15I2U|I are very similar to that appeared in mesonic sector Eq. 0. Therefore, we 
can employ the same method to evaluate these loop integrations. 
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FIG. 5: The scalar diquark-quark-quark coupling as a function of r^. The dotted and dash-dotted lines indicate the quark-quark 
pseudo-threshold for set A and B, respectively. 
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FIG. 6: The ratio of the on-shell approximation compared to the exact diquark T-matrix for the various scalar and axial 
diquark masses. 



We use the parameter sets determined in the mesonic sector shown in tabled Our numerical computation is valid 
below the first quark-quark pseudo-threshold. Note that the longitudinal polarizability Jail) defined in Eq. H20() does 
not vanish here. This term will be neglected in our one-pole approximation since it does not produce any poles in 
the T-matrix, and so makes a very small contribution compared to the transverse piece (see Fig.|2}. The longitudinal 
polarizability is not important in the local N JL model as well . We find that for a wide range of and , 

for all parameter sets, a bound scalar and axial-vector diquark exist (the results for additional sets can be found in 
|38j). This is in contrast to the normal NJL model where a bound axial- vector diquark exists only for very strong 
couplings The diquark masses for various values of and are plotted in Fig. 01 As already pointed out, the 
scalar diquark mass is equal to the pion mass at = 1. It is obvious from Fig. 21 that for — the axial- vector 
diquark is heavier than the scalar diquark, and consequently is rather loosely bound. For very small and r^, one 
finds no bound state in either diquark channels. In Fig. |SJ we show the scalar diquark-quark-quark coupling defined 
in Eq. (|17|l with respect to various scalar diquark couplings. 

One should note that the nucleon bound state in the diquark-quark picture does not require asymptotic-diquark 
states since the diquark state is merely an intermediate device which simplifies the three-body problem. Nevertheless, 
evidence for correlated diquark states in baryons is found in deep-inelastic lepton scatterings and in hyperon weak 
decays |39| . At the same time, diquarks appear as bound states in many phcnomcnological models, and are seen in 
lattice calculations Iii contrast to our perception of QCD colour confinement, the corresponding spectral 

functions for these supposedly confined objects in the colour anti-triplet channel are very similar to mesonic ones |4lj |. 

Next we study the off-shell behaviour of the diquark T-matrix. In Fig. El we show the discrepancy between the 
exact T-matrix and the on-shell approximation V^''^{q). At the pole we have by definition that V^'°'{q)\q2^j^,j2 = 1. 
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FIG. 7: The coupled Bethe-Salpeter equation for the vertex function 0. 



We see elsewhere that the off-shell contribution is very important due to the non-locality of our model. We find that 
the bigger the diquark mass is, the bigger the off-shell contribution. The off-shell behaviour of the scalar and the 
axial-vector channel for both parameter sets A and B are rather similar. 



IV. THREE-BODY SECTOR 



In order to make the three-body problem tractable, we discard any three-particle irreducible graphs. The relativistic 
Faddeev equation can be then written as an effective two-body BS equation for a quark and a diquark due to the 
locality of the form factor in momentum space (see Eq. and accordingly the separabilit^f the two-body interaction 
in momentum-space. Wc adopt the formulation developed by the Tiibingen group ,29, S^, Isij to solve the resulting 
BS equation. In the following we work in momentum space with Euclidean metric. The BS wave function for the 
octet baryons can be presented in terms of scalar and axialvector diquarks correlations, 

^(^'^H^'^) = (jI|5;;p))-(^'^)' (23) 

where u{P, s) is a basis of positive-energy Dirac spinors of spin s in the rest frame. The parameters p = {1 — ri)pi — 
ri{pj + pk) and P = pi + pj + pk are the relative and total momenta in the quark-diquark pair, respectively. The 
Mandelstam parameter rj describes how the total momentum of the nucleon P is distributed between quark and 
diquark. 

One may alternatively define the vertex function associated with ip{p, P) by amputating the external quark and 
diquark propagators (the legs) from the wave function; 

HP,P)-S-\p,)D-\p,)(^p^^''J^y (24) 

with 

where D{p), D^'^ (p) and S{p) are Euclidean versions of the diquark and quark propagators which are obtained by 
the standard transcription rules from the expressions in Minkowski space, Eqs. Hlt)l21|l and Eq. jS}, respectively. The 
spectator quark momentum pq and the diquark momentum pd are given by 

Pq = 7]P+p, (26) 
Pd = {l^v)P-P, (27) 

with similar expressions for kq^d, where we replace p by /c on the right-hand side. The Mandelstam parameter 77, 
parametrises different definitions of the relative momentum within the quark-diquark system. In the ladder approx- 
imation, the coupled system of BS equations for octet baryon wave functions and their vertex functions takes the 
compact form, 

p)-J -^.K^'^p^ ^; P)^ (28) 
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where K^^{p, k; P) denotes the kernel of the nucleon BS equation representing the exchange quark within the diquark 
with the spectator quark (see Fig.O, and in the colour singlet and isospin i channel we find (see Ref. ,26j) 

^ ' ' \-V5xHPi,kd)S^iQ)r{P2,Pd) -x"{Pi.kd)S^{<l)r{P2.Pd) )' ^ ' 

where x and (and their adjoint x ^^nd x'^) stand for the Dirac structures of the scalar and the axial- vector diquark- 
quark-quark vertices and can be read off immediately from Eqs. (|13l I16|l and Eqs. (|18l 121(1 . respectively. Therefore 
we have 



xHPi,kd) = gds,,h'C)^2V^{kd)f{pi + (1 - a)kd)f{-pi + akd), 

X^(pi,M = 9daqq{rC)^2V-{kd)f{pi + (1 - <y)kd)f{-pi + akd), (30) 

where cr is the Mandelstam parameter parameterising different definitions of the relative momentum within the 
quark-quark system. We have used an improved on-shell approximation for the contribution of diquark T-matrix 
occurring in the Faddccv equations. Instead of the exact diquark T-matrices we use the on-shell approximation with 
a correction of their off-shell contribution through V'''°'{p). What is neglected is then the contribution to the T- 
matrix beyond the pseudo-threshold. We will see this approximation is sufficient to obtain a three-body bound state. 
A similar approximation has been already employed in the normal NJL model in the nucleon sector, see for examples 
Refs. C)ne should note that here we do not have continuum states like the normal NJL model. However, 

there exist many complex poles beyond the pseudo-threshold which may be ignored, provided that they lie well above 
the energies of interest and the cutoff. For the parameter sets considered here, the next set of poles would result in 
another pseudo-threshold at energies of 1.3 GeV and 3 GeV for sets A and B, respectively. The model is not intended 
to be reliable at such momenta. On the other hand, as we will see in the next section, in practice, one may escape 
these poles far enough away by taking advantage of the above Mandelstam parametrization of the momenta. 

The relative momentum of quarks in the diquark vertices x a-nd sue defined as pi = p + fc/2 — (1 — 3ri)P/2 and 
P2 = —k — p/2 -|- (1 — 3?7)P/2, respectively. The momentum kd of the incoming diquark and the momentum pd of the 
outgoing diquark are defined in Eq. H27I) (see Fig. [T)). The momentum of the exchanged quark is fixed by momentum 
conservation ai q = —p — fc -I- (1 — 2rf)P. In the expressions for the momenta we have introduced two independent 
Mandelstam parameters 77, a, which can take any value in [0, 1]. Observables should not depend on these parameters 
if the formulation is Lorentz covariant. This means that for every BS solution ^{p, P] rji, ai) there exists an equivalent 
family of solutions. This provides a stringent check on calculations; see the next section for details. 

It is interesting to note that the non-locality of the diquark-quark-quark vertices naturally provides a regularisation 
of the ultraviolet divergence in the diquark-quark loop. 

We now constrain the Faddeev amplitude to describe a state of positive energy, positive parity and spin s = 1/2. 
The parity condition can be immediately reduced to a condition for the BS wave function: 



^ r{P,P) I 74A^V^(p,P)7V \-riP,p)'' ^ ' 



where we define p = K-pp and P — A-pP, with Aij^ = diag(— 1, —1, —1, 1). In order to ensure the positive energy 
condition, we project the BS wave function with the positive-energy projector A"*" = (1 -I- j^), where the hat denotes 
a unit four vector (in rest frame we have P = P/iM). Now we expand the BS wave function 'ijj(p, P) in Dirac 
space r g {1, 75, 7^, 757'', cr^''}. The above-mentioned conditions reduce the number of independent component from 
sixteen to eight, two for the scalar diquark channel. Si, (i = 1, 2) and six for the axial-diquark channel, Ai, (i — 1, ...6). 
The most general form of the BS wave function is given by 



ViP, P) = [iP'^T^tAi + P''A2 - py^As + ^pi^Ai + (pt^^j, ~ 7^) As 

-{ij^^t + Vt)A6)i5A+. (32) 



Here we write 7j^ = 7^ ^ /"P^. The subscript T denotes the component of a four- vector transverse to the nucleon 
momentum, Pt — P ^ P{p-P)- In the same way, one can expand the vertex function cf) in Dirac space, and since 
the same constraints apply to the vertex function, we obtain an expansion similar to Eq. (|32|l . with new unknown 
coefficients Si and A;. The unknown functions Si{^i) and Ai{Ki) depend on the two scalars which can be built from 
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the nucleon momentum P and relative momentum p, z = P.p = cosw (the cosine of the four-dimensional azimuthal 
angle of p^) and p^. Of course, they depend on as well, but this dependence becomes trivial in the nucleon rest 
frame. 

In the nucleon rest frame, one can rewrite the Faddecv amplitude in terms of tri-spinors each possessing definite 
orbital angular momentum and spin *2!^. It turns out that these tri-spinors can be written as linear combinations of 
the eight components defined in Eq. H^^2(l . Thus from knowledge of Si and Ai, a full partial wave decomposition can 
be immediately obtained [2^. Note that the off-shell contribution V'^'°-{q) is a function of the scalar q^. Moreover, 
the form factor in our model Lagrangian is also scalar, hence the total momentum dependent part of the diquark- 
quark-quark vertices are scalar functions and carry no orbital angular momentum, i.e. L'^x^'^il) — 0- Therefore, the 
partial wave decomposition obtained in Ref. p9l | for pointlike diquarks can be used here. Notice that no such partial 
wave decomposition can be found for the BS vertex function (f)^'^ since the axial-vector diquark propagator mixes the 
space component of the vertex function and time component of the axial-vector diquark. 



A. Numerical method for the coupled BS equations 

To solve the BS equations we use the algorithm introduced by Oettel et al l33 . The efficiency of this algorithm 
has already been reported in several publications, see for example Refs. [2^ ISOl iBlj . We will focus here only on the 
key ingredients of this method. The momentum dependence of quark mass in our model increases the complexity of 
the computation significantly. As usual, we work in the rest frame of the nucleon P = (0, iAIn)- In this frame we are 
free to chose the spatial part of the relative momentum p parallel to the third axis. Thus the momenta p and k are 
given by 



fc'" = |A:|(sin6l'sin(/)Vl - z'^, sin 61' cos - z'^,cos9\/l - z'^,z'), (33) 

where we write z = cosw and z' = cosw'. The wave function Eq. H32(l consists of 2 x 2 blocks in Dirac space can be 
simplified to 




Si{p\ z) \ , / a.VT^A.ip^z) 

asVT^S^ip^z) r ^ ' [ A2{p\z) 




i^/T^Ai{p^,z) Oy ' V \P. ) y -aiVT^Ae{p^,z) 

iaiAr,{p ^,z) o\ 
a^VT^Aeip^z) i ' 



(34) 



The great advantage of this representation is that the scalar and the axial- vector components are decoupled. Therefore 
the BS equation decomposes into two sets of coupled equations, two for the scalar diquark channel and six for the 
axial diquark channel. We expand the vertex (wave) functions in terms of Chebyshcv polynomials of the first kind, 
which are closely related to the expansion into hyperspherical harmonics. This decomposition turns out to be very 
efficient for such problems |23, Hfl Hll, Hj] . Explicitly, 

Ffip^z) = ^ z"^^f(")(p2)r„(z), 

Ff{p\z) = ^''Ff'-\p^)T^{z), (35) 

where Tn{z) is the Chebyshev polynomial of the first kind. We use a generic notation where the functions Ff{a.nd 
Pf) stand for any of the functions 5*^, Ai (and Si, A^), 

5*1,2 Pfa^ ^i...6 F3...S, 

§1,2 ^ fI^, A1...6 ^ fI„s. (36) 

We truncate the Chebyshev expansions involved in Pf and Pf at different orders nmax and mmax, respectively. We 
also expand the quark and diquark propagators into Chebyshev polynomials. In this way one can separate the P ■ p 
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and P ■ k dependence in Eqs. H24I28() . Using the orthogonality relation between the Chebyshev polynomials, one can 
reduce the four dimensional integral equation into a system of coupled one-dimensional equations. Therefore one can 
rewrite Eqs. H24I28|I in the matrix form 

j — 1 ni—O 
8 n„ 



' '-max pOO 

7 = 1 n=0 



(37) 



Here g"J^ and i?™" are the matrix elements of the propagator and the quark-exchange matrices, respectively. The 
indices n,m refer to the Chebyshev moments and i,j denote the individual channels. To solve Eq. (|37|l . we first 
rewrite it in the form of linear eigenvalue problem. Schematically 

X{P^)ip = K{P^)ip, (38) 

with the constraint that A(p2) = 1 at ^ -Mj^. This can be used to determine the nucleon mass Mn iteratively. 
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FIG. 8; The dependence of the nucleon mass on the Mandelstam parameter 77 for three choices of the cut-off on the Chebyshev 
expansion. Here we use set B, with Mds ~ 725 MeV and Alda ~ 630 MeV. The vertical lines at -qi and 7^2 denote the position 
of the singularities defined in Eq. Ijiffljl . 

As already pointed out, the BS solution should be independent of the Mandelstam parameters 77, a. As can be seen 
in Fig. 13 there is indeed a large plateau for the 77 dependence if we use a high cut-off on the Chebyshev moments. 
The limitations on the size of this area of stability can be understood by considering where the calculation contains 
singularities due to quark and diquark poles, 

,1 



7] G 



1 



1 - 



Mds 
Mn ' Mn 
Mda mi 



if 



if 



Mds < Mda, 



Mda < M, 



ds • 



(39) 



. The singularities in the quark-exchange propagator 
^^). No such constraint exists for a, which relates 



Mn Mn 

A similar plateau has been found in other applications |29ll30ll31l 

put another constraint on the acceptable range oi ij; rj > ^(1 — 
to the relative momentum between two quarks. To simplify the algebra we take a — 1/2. 

In what follows we use a momentum mesh of 60 x 60 for p, fc, mapped in a non-linear way to a finite interval. In the 
non-singular regime of Mandelstam parameter 77 Eq. (|39|l , the Faddeev solution is almost independent of the upper 
limit on the Chebyshev expansion, and for m,nax — 10,7imax — 12, see the Fig.lHl this seems to be satisfied. This limit 
is somewhat higher than the reported values for simple models |2 



B. Nucleon Solution 



In order to understand the role of the axial diquark in nucleon solution, we first consider the choice = 0. For 
this case we find that the non-confining set A can not generate a three-body bound state without the inclusion of the 
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FIG. 9: The nucleon mass without inclusion of the axial diquark channel. The dotted lines indicate the diquark-quark threshold. 
All values are given in GeV. 
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FIG. 10: The nucleon mass as a function of the scalar diquark mass for various axial vector diquark masses for both parameter 
sets. The scalar diquark-quark threshold are shown by the dotted lines. 



ofF-shell contribution. For the confining set B one also has to enhance the diquark-quark-quark couphng gasqq by a 
factor of about 1.73 over the value defined in Eq. p7|l (as we will show, this extra factor is not necessary when the 
axial- vector diquark is included) . The situation is even more severe in the on-shell treatment of the local N JL model, 
since one needs to include the quark-quark continuum contribution in order to find a three-body bound state when 
the axial- vector diquark channel is not taken into account ^3|- 

As can be seen from Fig. 0] a decrease in rs leads to a larger diquark mass, and an increase in the off-shell 
contribution to the quark-quark T-matrix (see Fig.|51). This off-shell correction is crucial for forming a bound nucleon. 
If the off-shell contribution V^{q) is omitted, a bound nucleon can not be found. 

The nucleon result is shown in Fig.|51 We also show a fictitious diquark-quark threshold defined as Mds + ti^. The 
nucleon mass can be seen to depend roughly linearly on the scalar diquark mass. A similar behaviour is also seen in 
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FIG. 11: Range of parameters {rs,ra) where we find a nucleon mass of 940 MeV. 



TABLE III: Diquark masses and coupling of diquarks to quarks obtained for Mm = 940 MeV. All masses are given in MeV. 
Eda{Eda) denote the binding energy of diquarks in the nucleon, Eff{Ef^) denote the binding energy of the nucleon measured 
from scalar (axial) diquark mass. 
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Set A2 


Set A3 


Set Bl 


Set B2 


Set B3 


Mds 


775 


748 


698 


802 


705 


609 


Qdsqq 


0.74 


0.83 


1.04 


0.73 


1.31 


1.79 


rs 


0.09 


0.12 


0.17 


0.06 


0.14 


0.24 


Eda 


7 


34 


84 


14 


111 


207 




226 


199 


149 


270 


173 


77 


Mda 


705 


725 


775 


604 


660 


725 


Qdaqq 


1.08 


0.98 


0.79 


1.99 


1.67 


1.28 


Ta 


0.20 


0.17 


0.11 


0.32 


0.23 


0.15 


Eda 


77 


57 


7 


212 


156 


91 


T-^da 


156 


176 


226 


72 


123 


193 


RMS 


194.88 


181.51 


163.70 


283.99 


232.86 


209.95 



the local NJL model j23. Increasing the diquark mass (or decreasing rs) increases the nucleon mass, i.e. the scalar 
diquark channel is attractive. In order to obtain a nucleon mass of 940 MeV, we need diquark mases of 608 MeV and 
623 MeV for set A and B, respectively. The corresponding nucleon binding energy measured from the diquark-quark 
threshold are 56 MeV and 91 MeV for set A and B, respectively, compared to the binding of the diquarks (relative to 
the quark-quark pseudo threshold) of about 174 and 193 MeV for set A and B, respectively. Such diquark clustering 
within the nucleon is also observed in the local NJL model [25| . and is qualitatively in agreement with an instanton 
model [i^ and lattice simulations |4^ . 

Next we investigate the effect of the axial- vector diquark channel on nucleon solution. We find that the axial- vector 
diquark channel contributes considerably to the nucleon mass and takes away the need for the artificial enhancement 
of the coupling strength for set B. In Fig. ^|we show the nucleon mass as a function of the scalar and axial- vector 
diquark mass. As in the scalar diquark channel, we define the axial-vector diquark-quark threshold as Mda + 
We see that as one increases the axial-vector diquark (and scalar diquark) masses, the quark-quark interaction is 
weakened and consequently the nucleon mass increases. Therefore the contribution of the axial-vector channel to the 
nucleon mass is also attractive. 

In Fig. 111! we plot the parameter space of the interaction Lagrangian with variable r^. and which leads to the 
nucleon mass = 940 MeV. The trend of this plot for the non-confining set A is very similar to the one obtained 
in the local NJL model 26] . 

If the scalar diquark interaction is less than 0.14, we need the axial- vector interaction to be stronger than the 
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FIG. 12: Chebysliev moments (labelled by n) of the scalar and axial-vector (AV) diquark amplitudes of the nucleon BS wave 
function given by Set Al. 
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FIG. 13: diebyshev moments (labelled by the order n) of scalar and axial- vector (AV) diquark amplitudes for the nucleon BS 
wave function obtained for parameter set Bl. 
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FIG. 14: Diagram corresponding to the definition of density Eq. H4()|l . 

scalar diquark channel > rs in order to get the experimental value of nucleon mass. For set B, as we approach 
ra = 0, the curve bends upward, reflecting the fact that we have no bound state with only the scalar diquark channel. 

In Fig. ^2 we see for the confining set B that the interaction is again shared between the scalar and the axial- vector 
diquark and for small < 0.19 one needs a dominant axial- vector diquark channel ra > Vg- It is obvious that the 
axial-vector diquark channel is much more important in the confining than the non-confining phase of model. 

In order to study the implications of the quark confinement for the description of the nucleon, we compare in Table 
mil three representative cases for both the non-confining and confining parameter sets, which all give a nucleon mass 
of about 940 MeV. The first three columns contain results for set and the last three columns for the confining set 
B. 

Given the definition of diquark-quark thresholds, in the presence of both scalar and axial-vector diquark channels, 
the diquarks in the nucleon can be found much more loosely bound, although one obtains a very strongly bound 
nucleon solution near its experimental value, see table lllll Next we study the nucleon BS wave function for the 
various sets given in Table Uni The nucleon wave and vertex function are not physical observables, but rather they 
suggest how observables in this model will behave. In Figs. El El we show the leading Chebyshev moments of the 
scalar functions of the nucleon BS wave function for various sets (Al and Bl). They describe the strengths of the 
quark-diquark partial waves with S" as a total quark-diquark spin and L as a total orbital angular momentum. They 
are normalised to Ff'''^^{pi) — 1, where pi is the first point of the momentum mesh. Very similar plots are found for 
the other sets of parameters given in table IIIII It is seen that the contribution of higher moments are considerably 
smaller than lower ones, indicating a rapid convergence of the expansion in terms of Chebyshev polynomials. In the 
confining case. Fig. 1131 there is a clear interference which is not present in the non-confining one. Fig. 1121 Therefore 
in the confining case, all wave function amplitudes are shifted to higher relative four-momenta between the diquark 
and quark. 

In order to understand the effect of this interference, we construct a density function for the various channels in 
the nucleon rest frame. This density is defined as 

piP±,P)^ J dpi^\p^,Pi,P)D-\pdMp^,P4,P), (40) 

where p± stands for the space component of the relative momentum p, and D^^{pd) is defined in Eq. I|25|l . This 
definition corresponds to a very naive diagram describing the quark density within the nucleon, see Fig. 1141 In the 
above definition of the density function, we have integrated over the time component of the relative momentum. 
In this way the density function becomes very similar to its counterpart in Minkowski space. Although the above 
definition of density is not unique, it does provide a useful measure of the spatial extent of the wave function. Similar 
calculations have been done for the quark condensate in Ref. [23|. The results are plotted in Figs. 1151 andllGI 

It is noticeable that in all cases, the s-wave in scalar diquark channel is the dominant contribution to the ground 
state. The relative importance of the scalar and the axial diquark amplitude in the nucleon changes with the strength 
of the diquark-quark couplings gdsqq{gdaqq) and accordingly with rs{ra)- There are indications that in the confining 
sets, the nucleon density extends to higher relative momentum between the diquark and the quark. These imply a 
more compact nucleon in the confining cases. In order to find a quantitative estimate of the confinement effect in our 
model, we calculate p^^^ — [{p\) — {p±)^y^^; the results can be found in table IIIll 

We also see in the both confining and non-confining cases a decrease in pf^^'^s .^[f^iy weakening axial-vector diquark 
interaction (and consequently increasing the scalar diquark interaction strength). This can be associated with the 
important role of the axial-vector diquarks. If we compare p^^ for the two sets A2 and B2, which have very similar 
interaction parameters rs(ra), an increase about 25% is found. 
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FIG. 15: Shows the nucleon density [M = 940 MeV, Set A) with respect to relative momentum between diquark and quark 
for different set of A1,A2 and A?,. 



V. SUMMARY AND OUTLOOK 



In this work we have investigated the two- and three-quark problems in a non-local NJL model. We have truncated 
the diquark sector to the scalar and the axial-vector channels. We have solved the relativistic Faddeev equation for 
this model and have studied the behaviour of the nucleon solutions with respect to various scalar and the axial- vector 
interactions. We have studied the dependence of the baryon masses and waves on the interaction parameters and 
(and on the scalar and the axial- vector diquark mases), which can describe the pion and the nucleon simultaneously. 
In order to put more restrictions on these parameters, one would need to calculate the A mass as well. 

Although the model is quark confining, it is not diquark confining (at least in the rainbow-ladder approximation). 
A bound diquark can be found in both scalar and the axial-vector channel for a wide range of couplings. We have 
found that the off-shell contribution to the diquark T-matrix is crucial for the calculation of the nucleon: without 
it the attraction in the diquark channels is too weak to form a three-body bound state. We have also found that 
both the scalar and the axial- vector contribute attractively to the nucleon mass. The role of axial- vector channel is 
much more important in the confining phase of model. The nucleon in this model is strongly bound even though the 
diquarks are rather loosely bound. The confining aspects of the model are more obvious in three-body, rather than 
the two-body sector. We decomposed the nucleon BS wave function in the nucleon rest frame in terms of spin and 
orbital angular momentum eigenstates and constructed the quark density function in the various channels. It was 
revealed that s-wave in scalar diquark channel is the dominant contribution to the ground state. By investigating 
the nucleon wave function we found that quark confinement leads to a more compact nucleon. The size of nucleon is 
reduced by about 25% in the confining cases. 
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FIG. 16: Shows the nucleon density (M = 940 MeV, Set B) with respect to relative momentum between diquark and quark 
for different set of Bl, B2 and B3. 



For both confining and non-confining cases, an increase in the scalar diquark channel interaction leads to a 
lower nucleon mass , se e Fig. 1101 However, the mass of the A should be independent of since it does not contain 
scalar diquarks [2^ [23, HI] . In the standard NJL model the difference between the nucleon and A masses is strongly 
dependent on the scalar diquark interaction . In the current model where the axial- vector diquark makes a larger 
contribution to the nucleon mass, a detailed calculation of the A is needed to understand the mechanism behind the 
A-7V mass difference. Note also that neglecting TriV-loops may lead to a quantitative overestimate of the axial- vector 
diquark role in the nucleon f45l |. 

In order to understand the implications of this model in baryonic sector fully one should investigate other properties 
of the nucleons such as the charge radii, magnetic moments and axial coupling. On the other hand, the role of quark 
confinement in this model may be better clarified by investigating quark and nuclear matter in this model. Such 
problems can also be studied within the same Faddeev approach 0, |4g . 
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